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Quantum memories for light, which allow the reversible transfer of quantum states between light 
and matter, are central to the development of quantum repeaters pQ, quantum networks [2], and 
linear optics quantum computing [3]. Significant progress has been reported in recent years, includ- 
ing the faithful transfer of quantum information from photons in pure and entangled qubit states 
[4l410j . However, none of these demonstrations confirm that photons stored in and recalled from 
quantum memories remain suitable for two-photon interference measurements, such as C-NOT gates 
and Bell-state measurements, which constitute another key ingredient for all aforementioned appli- 
cations of quantum information processing. Using pairs of weak laser pulses, each containing less 
than one photon on average, we demonstrate two-photon interference as well as a Bell-state measure- 
ment after either none, one, or both pulses have been reversibly mapped to separate thulium-doped 
titanium-indiffused lithium niobate (Ti:Tm:LiNb03) waveguides. As the interference is always near 
the theoretical maximum, we conclude that our solid-state quantum memories, in addition to faith- 
fully mapping quantum information, also preserves the entire photonic wavefunction. Hence, we 
demonstrate that our memories are generally suitable for use in advanced applications of quantum 
information processing that require two-photon interference. 



When two indistinguishable single photons impinge on 
a 50/50 beam-splitter (BS) from different input ports, 
they bunch and leave together by the same output port. 
This so-called Hong-Ou-Mandel (HOM) effect pi] is due 
to destructive interference between the probability am- 
plitudes associated with both input photons being trans- 
mitted or both reflected, see Fig. 1. Since no such in- 
terference occurs for distinguishable input photons, the 
interference visibility V provides a convenient way to ver- 
ify that two photons are indistinguishable in all degrees 
of freedom, i.e. spatial, temporal, spectral, and polariza- 
tion modes. The visibility is defined as 

^ (^max ^-min)/'^'max5 (1) 

where 7?. m i n and 7?. max denote the rate with which pho- 
tons are detected in the two output ports in coincidence 
if the incoming photons are indistinguishable and distin- 
guishable, respectively. Consequently, the HOM effect 
has been employed to characterize the indistinguishabil- 
ity of photons emitted from a variety of sources, including 
parametric down-conversion crystals [12] . trapped neu- 
tral atoms [T3] , trapped ions [2] , quantum dots [T51 - 
IT7] . organic molecules [18], nitrogen- vacancy centres in 
diamond [19j [20] , and atomic vapours [2TH25] . Further- 
more, two-photon interference is at the heart of linear 
optics Bell-state measurements [55] , and, as such, has al- 
ready enabled experimental quantum dense coding [2"T] . 
quantum teleportation [55], and entanglement swapping 
[2"5] , However, to date, the possibility to perform Bell- 
state measurements with photons that have previously 
been stored in a quantum memory, as required for ad- 
vanced applications of quantum information processing, 



has not yet been established. For these measurements 
to succeed, photons need to remain indistinguishable in 
all degrees of freedom, which is more restrictive than 
the faithful recall of encoded quantum information. In- 
deed, taking into account that photons may or may not 
have been stored before the measurement, this criterion 
amounts to the requirement that a quantum memory pre- 
serves a photon's wavefunction during storage. Similar to 
the case of photon sources, the criterion of indistinguisha- 
bility is best assessed using HOM interference, provided 
single-photon detectors are employed. 

Our experimental setup is depicted in Fig. [2] We 
employ solid-state quantum memories, more precisely 
thulium-doped lithium-niobate waveguides in conjunc- 
tion with the atomic frequency comb (AFC) quantum 
memory protocol [30] . which have shown great promise 
for advanced applications of quantum information pro- 
cessing [5HZ]- We then interfere various combinations of 
recalled and non-stored (i.e. directly transmitted) pulses 
on a 50/50 BS (HOM-BS). When using single photon 
Fock states at the memory inputs, the HOM visibility 
given in Eq. theoretically reaches 100% as illustrated 
in Fig. [l] However, with phase incoherent laser pulses 
obeying Poissonian photon-number statistics, as in our 
demonstration, the maximally achievable visibility is 50% 
[31) . irrespective of the mean photon number (see Sup- 
plementary Information). Nevertheless, attenuated laser 
pulses are perfectly suitable for assessing the effect of our 
quantum memories on the photonic wavefunction. Any 
reduction of indistinguishability due to storage causes a 
reduction of visibility, albeit from maximally 50%. This 
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FIG. 1: Illustration of HOM-interference in the case of single photons at BS input \tl> ( - ln ^ 1 ) = j 1, 1'}, where the prime on the latter 
input indicates the possibility to distinguish that input photon from the other in some degree of freedom e.g. by being polarized 
orthogonally. The four possible paths of the photons are illustrated, together with their corresponding output states. If the 
input photons are indistinguishable with respect to all degrees of freedom we can ignore the primes in the output states and 
the paths shown in the two central pictures are identical and, due to the different signs, thus cancel. This leaves in the output 
state |7/^ out '} only the possibilities in which photons bunch. For distinguishable photons, e.g. having orthogonal polarizations, 
all paths are distinguishable and all terms remain in |t/> ( ' out - ) }. 



approach extends the characterization of quantum mem- 
ories using attenuated laser pulses [35] from assessing the 
preservation of quantum information during storage to 
assessing the preservation of the entire wavefunction, and 
from first- to second-order interference. 

We first deactivate both quantum memories (see Sup- 
plementary Information), to examine the interference be- 
tween directly transmitted pulses, and thereby establish 
a reference visibility for our experimental setup. We set 
the mean photon number per pulse before the memories 
to 0.6, i.e. to the single-photon level. Using the wave 
plates, we rotate the polarizations of the pulses at the 
two HOM-BS inputs to be parallel (indistinguishable) or 
orthogonal (distinguishable). Employing Eq. ([I]) we find 
a visibility of (47.9 ±3.1)%. 

Subsequently, we activate memory a while keeping 
memory b off, and adjust the timing of the pulse prepa- 
ration so as to interfere a recalled pulse from the active 
memory with a directly transmitted pulse from the in- 
active memory (see Supplementary Information). Pulses 
are stored for 30 ns in memory a, and the mean photon 
number per pulse at the quantum memory input is 0.6. 
Taking the limited storage efficiency of « 1.5% and cou- 
pling loss into account, this results in 3.4 x 10 -4 photons 
per pulse at the HOM-BS inputs. As before, changing the 
pulse polarizations from mutually parallel to orthogonal, 
we find V = (47.7 ± 5.4)%, which equals our reference 
value within the measurement uncertainties. 

As the final step, we activate both memories to test 
the feasibility of two-photon interference in a quantum- 
repeater scenario. We note that in a real-world im- 
plementation, memories belonging to different network 



nodes are not necessarily identical in terms of material 
properties and environment. This is captured by our 
setup where the two Ti:Tm:LiNb03 waveguides feature 
different optical depths and experience different magnetic 
fields (see Fig. [2] and Supplementary Information) . To 
balance the ensuing difference in memory efficiency we 
set the mean photon number per pulse before the less 
efficient and more efficient memories to 4.6 and 0.6, re- 
spectively, so that, as before, the mean photon num- 
bers are 3.4 x 10~ 4 at both HOM-BS inputs. With 
the storage time of both memories set to 30 ns, we get 
V = (47.2 ± 3.4)%, in excellent agreement with the val- 
ues from the previous measurements. The consistently 
high visibilities, compiled in the first column of Table [TJ 
hence confirm that our storage devices do not introduce 
any degradation of photon indistinguishability during the 
reversible mapping process, and that two-photon inter- 
ference is feasible with photons recalled from separate 
quantum memories, even if the memories are different. 

We now investigate in greater detail the change in coin- 
cidence count rates as photons gradually change from be- 
ing mutually indistinguishable to completely distinguish- 
able w.r.t. each degree of freedom accessible for change 
in single-mode fibres, i.e. polarization, temporal, and 
spectral modes (see Supplementary Information). To ac- 
quire data more efficiently we increase the mean number 
of photons per pulse at the memory input to between 10 
and 50 (referred to as few-photon-level measurements). 
However, the mean photon number at the HOM-BS re- 
mains below one. Example data plots are shown in Fig.[3j 
while the complete set of plots is supplied in the Supple- 
mentary Information Figs. 6][£ 
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FIG. 2: Experimental setup. Light from a 795.43 nm wavelength CW laser passes through an acousto-optic modulator 
(AOM) driven by a sinusoidally varying signal. The first negative refraction order is fibre coupled into a phase modulator 
and, via a beam-splitter (BS), two polarization controllers (PCs) and two micro-electromechanical switches (MEMS), injected 
from the back into two Ti:Tm:LiNb03 waveguides (labelled a and b) cooled to 3 K 33 . Waveguide a is placed inside a 
superconducting solenoid. Using a linear frequency-chirping technique [34] we tailor AFCs with 600 MHz bandwidth and a 
few tens of MHz peak spacing, depending on the experiment, into the inhomogeneously broadened absorption spectrum of 
the thulium ions, as shown for crystal a in the inset. After 3 ms memory preparation time and 2 ms wait time we store and 
recall probe pulses during 3 ms. The 8 ns long probe pulses with « 50 MHz Fourier-limited bandwidth are derived from the 
first positive diffraction order of the AOM output at a repetition rate of 2.5-3 MHz. Each pulse is divided into two spatial 
modes by a half- wave plate (HWP) followed by a polarizing beam-splitter (PBS). All pulses are attenuated by neutral-density 
filters (NDFs) and coupled into optical fibres and injected from the front into the Ti:Tm:LiNb03 waveguides. After exiting the 
memories (i.e. either after storage, or after transmission), the pulses pass quarter- and half- wave plates used to control their 
polarizations at the 50/50 BS (HOM-BS) where the two-photon interference occurs. Note that, to avoid first-order interference, 
pulses passing through memory a propagate through a 10 km fibre to delay them w.r.t. the pulses passing through memory b 
by more than the laser coherence length. Finally, they are detected by two single-photon detectors (actively quenched silicon 
avalanche photodiodes, Si-APDs) placed at the outputs of the beam-splitter, and coincidence detection events are analyzed 
with a time-to-digital convertor (TDC) and a computer. 



In Fig. [3ji we show the coincidence counts rates as a 
function of the polarization of the recalled pulse for the 
case of one active memory. The visibilities for all configu- 
rations (i.e. zero, one, or two active memories) extracted 
from fits to the experimental data are listed in column 2 
of Table [TJ They are - as in the case of single-photon-level 
inputs - equal to within the experimental uncertainty. 

Next, in Fig. [3]d, we depict the coincidence count rates 
as a function of the temporal overlap (adjusted by the 
timing of the pulse generation) for the two-memory con- 
figuration. Column 3 of Table |T] shows the visibilities 
extracted from Gaussian fits to the data, reflecting the 
temporal profiles of the probe pulses, for all configura- 
tions. Within experimental uncertainty, they are equal 
to each other. Alternatively, in the single-memory con- 
figuration, we also change the temporal mode overlap 
by adjusting the storage time of the pulse mapped to 
the quantum memory. Again the measured visibility of 
V = (44.4 ±6.9)% (see Fig.|3j:) is close to the theoretical 
maximum. 

Finally, we vary the frequency difference between the 
two pulses (see Supplementary Information) to witness 
two-photon interference w.r.t. spectral distinguishability. 
For this measurement, we consider only the configura- 
tions in which neither, or a just single memory is active. 
In both cases the visibilities, listed in the last column of 
Table [IJ are around 43%. While this is below the visibil- 



ities found previously, for reasons discussed in the Sup- 
plementary Information, the key observation is that the 
quantum memory does not affect the visibility. 



TABLE I: Experimental two-photon interference visibilities 
(%) for different degrees of freedom 





Single-photon 


Few-photon 


Storage 


level 


level 


configuration 


Polarization 


Polarization Temporal Spectral 


No-storage 


47.9 ±3.1 


51.0 ±5.6 42.4 ± 2.3 43.7 ± 1.7 


Single-storage 


47.7 ±5.4 


55.5 ± 4.1 47.6 ± 3.0 42.4 ± 3.5 


Double-storage 


47.2 ±3.4 


53.1 ±5.3 46.1 ±3.2 N. A. 



As stated in the introduction, Bell-state measurements 
(BSM) with photonic qubits recalled from separate quan- 
tum memories are key ingredients for advanced applica- 
tions of quantum communication. To demonstrate this 
important element, we consider the asymmetric (and ar- 
guably least favourable) case in which only one of the 
qubits is stored and recalled. Appropriately driving the 
AOM in Fig. [2] we prepare the states and 1^2)1 

which describe time-bin qubits [35] of the form |e), \l), 
^7j(l e ) + 10)' or (|e) — |Z)), where e and I, respectively, 
label photons in early or late temporal modes, which are 
separated by 25 ns. The qubits are directed to the mem- 
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FIG. 3: HOM interference plot examples for one or two active memory configurations (as labelled), a) Varying mutual 
polarization difference, b) Varying temporal overlap by changing timing of pulse generation, c) Varying temporal overlap by 
changing storage time. The acquisition time per data point is 60 s in a,b and 120 s in c. 



ories of which only one is activated. The mean photon 
number of the qubit that is stored is set to 0.6, yielding 
a mean photon number of both qubits at the HOM-BS 
input of 6.7 x 10 -4 . We ensure to overlap pulses encoding 
the states and |^) a t the HOM-BS and count coin- 
cidence detections that correspond to a projection onto 
the \i/>~) = ^(|e)|0 - |/)|e)) Bell state. This projection 
occurs if the two detectors click with 25 ns time difference 
[35j . Because \ip~) is antisymmetric w.r.t. any basis, the 
count rate is expected to reach a minimum value 1Z\\ if 
the two input pulses are prepared in equal states, and 
a maximum value 1Z± if prepared in orthogonal states. 
Accordingly, we define an error rate that quantifies the 
deviation of the minimum count rate from its ideal value 
of zero: 



ft" 



ft" +TZ 1 - 



(2) 



First, choosing to encode and jvE^) in states |e) and 



\l) we obtain the error rate — 0.039 ± 0.037, which 

is near the theoretical value of efj^ = (see the Sup- 
plementary Information for derivations of the theoretical 
values and bounds). In addition it clearly violates the 



lower bound e 



(CM) 

:/l 



0.33 that can be obtained for a 



Bell-state measurement on two qubits of which one is 
recalled from a classical memory (CM). Note that val- 
ues for and efjf 1 ^ are independent of whether |e) 
and \l) qubits are encoded into single photons or at- 
tenuated laser pulses. Next, using instead the states 



1+)^ 
c ( ex p) 
-+/- 



) + |^)), and |— ) = 4§(|e) — 10) we measure 
= 0.287 ± 0.020, which again only slightly exceeds 



75* 



the lowest possible value for attenuated laser pulses of 
e (att,QM) _ q 25 cruc j a i observation is once more 

that e^?_ violates both the lower bound for qubits en- 



coded into single photons e ^" l ^' CM ^ 



ated laser pulses e 



(att,CM) 

+ /- 



0.33 and attenu- 



0.417 that are recalled from 



a classical memory. 

Our demonstrations show that solid-state AFC quan- 
tum memories are suitable for two-photon interference 
experiments, even in the general case of storing the two 
photons an unequal number of times. With improved sys- 
tem efficiency |36j and multi-mode storage supplemented 
by read-out on demand |37H39) . such memories can be 
used as synchronization devices in multi-photon experi- 
ments, which will allow increasing the number of photons 
that can be harnessed simultaneously for quantum infor- 
mation processing beyond the current limit of eight [40j . 
A subsequent goal is to develop workable quantum re- 
peaters or, more generally, quantum networks, for which 
longer storage times are additionally needed. Depend- 
ing on the required value, which may range from hun- 
dred micro-seconds [41] to seconds [U [42] , this may be 
achieved by storing quantum information in optical co- 
herence, or it may require mapping of optical coherence 
onto spin states [50] . 
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SUPPLEMENTARY INFORMATION 
Memory operation and properties. 

A quantum memory is said to be activated when we 
configure the MEMS to allow the optical pumping light 
to reach the waveguide during the preparation stage and 
thus tailor an AFC in the inhomogeneously broadend ab- 
sorption spectrum of thulium ions (see Fig. [5]). If the 
optical pumping is blocked, the memory is said to be de- 
activated and light entering the waveguide merely experi- 
ences constant attenuation over its entire spectrum. If a 
memory is activated, an incident photon is mapped onto 
a collective excitation of thulium ions in the prepared 
AFC and subsequently re-emitted at a time given by the 
inverse of the comb tooth spacing [3U], i.e., t = 1/A (see 
Fig. [2]). In all cases, we adjust the mean photon number 
at the memory inputs so that mean photon numbers are 
equal at the HOM-BS inputs. This is required for achiev- 
ing maximum visibility with attenuated laser pulses (fur- 
ther details later in Supplementary Information). 

The two Ti:Tm:LiNb03 waveguides are fabricated 
identically but differ in terms of overall length, yielding 
optical depths of 2.5 for memory a and 3.2 for memory b. 
As shown in Figure 1, memory a is placed at the centre 
of a solenoid in a uniform magnetic field, while memory b 
is placed outside the solenoid and thus experiences only 
a much weaker stray field. Therefore it is not possible to 
achieve the optimal efficiency for both memories at the 
same time (further details later in Supplementary Infor- 
mation). 

Changing degrees of freedom. 

a) The polarization degree is easily adjusted using 
the free-space half- and quarter-wave plate set at each 
HOM-BS input. For our measurements we rotate the 
half-wave plate in steps of either 45 or 7.5 degrees, b) 
The temporal separation St between a pulse arriving at 
one of the HOM-BS inputs and the next pulse in the 
train arriving at the other input can be expressed as 
St = {?ii/c} mod St r , where n is the refractive index of 
the fibres, I ~ 10 km is the path-length difference for 
pulses interacting with memory a and 6, and St r is the 
repetition period of the pulse train from the AOM, which 
is set in the range of 350-400 ns. As we can change 5t r 
with 10 ps precision, we can tune St on the ns scale, c) 
For the storage time scan, the recall efficiency decreases 
with storage time due to decoherence. Hence, we balance 
the mean photon number per pulse for stored and trans- 
mitted pulses for each storage time, d) Finally, to change 
the spectral overlap of the pulses input to the HOM-BS 
we can utilize that these pulses were generated at dif- 
ferent times in the AOM and thus we can chose their 
carrier frequencies independently. We interchangeably 



drive the AOM by frequencies v a and and thus create 
two interlaced trains of pulses with different frequencies. 
Adjusting the pulse timing we can ensure that the pulses 
overlapped at the HOM-BS belong to different trains and 
thus have a spectral overlap given by Sv = v a — z/&. Due 
to the limited bandwidth of the AOM we are only able 
to scan Sv by 100 MHz, which, when compared to the 
50 MHz pulse bandwidth, is not quite sufficient to make 
the pulses completely distinguishable. To achieve com- 
plete distinguishability, we supplement with a measure- 
ment using orthogonal polarizations at the inputs (fur- 
ther details later in Supplementary Information). 



Preparing states for Bell-state measurement. 

For the Bell-state projection measurement we inter- 
changeably prepare the time-bin qubits in either |e) or 
|Z), or in i(|e) + \l)) and ^(|e) - \l)) by setting the 
relative phase and intensity of the AOM drive signal. 
Adjusting the timing of the pulse preparation we ensure 
that qubits in different states overlap at the HOM-BS. 



Properties of waveguide LiNb03 crystal and AFC. 

In the experimental configuration in which the HOM- 
interference occurs between two pulses recalled from sep- 
arate quantum memories we pointed, in the main text, 
to the different properties of the two memory devices. In 
this section we wish to elaborate on the differences be- 
tween the two memories based on their physical dissim- 
ilarity and measured optical depth as a function of fre- 
quency Memory waveguide a is 10.4 mm long and crystal 
b is 15.4 mm long. The optical depths at 795.43 nm are 
around 2.5 and 3.2 for waveguide a and b, respectively, as 
shown by the light-grey curves in Fig. [4] a, b, correspond- 
ing to the case in which the memories are not activated. 

In order to spectrally tailor an AFC in Tm:LiNb03, a 
magnetic field must be applied along the crystal's c-axis 
so as to split the ground and excited level multiplets into 
their two nuclear Zeeman sublevels [33j . However, as one 
crystal is located at the centre of the setup's solenoid 
and the other outside the solenoid (see Fig. [2]) it is not 
possible to apply the same B-field at the two crystals. 
Thus when we activate both memories we generally ap- 
ply a magnetic field, which provides a reasonable bal- 
ance in recall efficiencies but is not optimal for either 
memory. This circumstance is reflected by the different 
shapes of optical-depth profiles of the AFCs shown in red 
in Fig. [4ji,b. 
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FIG. 4: Measured optical depths of our two Ti:Tm:LiNb03 waveguides as a function of frequency shift of the probing light 
imparted by the phase-modulator. Light grey traces show optical depths when the memories are inactive, i.e. no AFC is 
prepared. Dark red traces show the prepared AFCs at a magnetic field of 900 Gauss at the centre of the solenoid. 



Two-photon interference in imperfectly prepared 
memories. 

In all our demonstrations of the HOM interference we 
consistently observe that the HOM visibility is close to 
the theoretical maximum for coherent states. Yet, it is 
important to realize that an improperly configured AFC 
quantum memory does alter a stored photon's wavefunc- 
tion, resulting in imperfect HOM interference with a non- 
stored photon. 

To support this claim we activate only memory a, 
whose performance we change by varying the bandwidth 
of the AFC, and interfere the recalled pulses with pulses 
directly transmitted through the deactivated memory b. 
As the AFC bandwidth decreases below that of the probe 
pulses, the AFC effectively acts as a bandpass filter for 
the stored photons and we thus expect the recalled pulses 
to be temporally broadened w.r.t. the original pulse. 
This is observed in the insert of Fig. [5j which shows 
smoothed histograms of photon detection events as a 
function time. It is worth noting that the small band- 
width AFC also acts as a bandpass filter for the trans- 
mitted pulse by virtue of the different effective optical 
depths inside and outside the AFC. Thus the broadened 
transmitted pulse starts to overlap with the echo for the 
narrow AFC bandwidth traces, as is also observed in the 
insert of Fig. [5] 

Another consequence of reducing the AFC bandwidth 
is that the overall efficiency of the quantum memory de- 
creases, which causes an imbalance between the mean 
photon numbers at the HOM-BS inputs and thus reduces 
HOM interference visibility. We circumvent the change 
to the echo efficiency by adapting the mean photon num- 



ber at the memory input so as to keep the mean photon 
number of the recalled pulse constant. With this reme- 
dial procedure, we assess the HOM visibility by changing 
the HOM-BS inputs from parallel to orthogonal polar- 
izations for a scries of different AFC bandwidths. The 
HOM visibility in Fig. [5] is steady for bandwidths from 
around 100 MHz and up. However, below 100 MHz the 
visibility begins to drop significantly. The dashed line is 
a fit of the visibilities to a Gaussian function with full- 
width at half-maximum (FWHM) of 79 ± 4 MHz. Note, 
that the reason for the visibility being limited to around 
40% is solely that, for this measurement, we do not go 
through the usual careful optimization steps. 

With these measurements we have illustrated how a 
quantum memory could alter the photonic wavefunction 
resulting in a reduced HOM interference visibility. A 
combination of spectral and temporal distortion of the 
photonic wavefunction is indeed a common type of per- 
turbation by quantum memories. FSJ |U] It is particu- 
larly worth noting that the gradient-echo memory (GEM) 
quantum memory protocol, though similar to the AFC 
protocol, imparts a frequency chirp to the recalled pulse 
[45] . If not corrected, this feature constitutes a pertur- 
bation of the wavefunction of the recalled pulse, which 
may render it unsuitable for applications relying on two- 
photon interference. 



Analytical model of second-order interference in 
coincidence measurements. 

In the following theoretical treatment we will derive ex- 
pressions for the coincidence and single-detector counts 
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FIG. 5: HOM interference visibility if HOM-BS input pulses are recalled from AFCs with varying bandwidths. Insert: His- 
tograms of recalled pulse detection times for different AFC bandwidths clearly showing broadening of recalled (and transmitted) 
pulses for bandwidths below 100 ns. 



in terms of probabilities. By multiplying these probabili- 
ties with the average experimental repetition rate we can 
easily calculate the predicted experimental count rates. 
To a large extent though, we will mainly be interested 
in relative probabilities or count rates between different 
settings of the degrees of freedom of pulses. 

It is reasonably straightforward to derive the rates of 
detection of photons at the outputs of a BS (note that 
in this Supplementary Information, the HOM-BS of the 
main text will be referred to as just BS) In our case co- 
herent states | a) and \f3), characterized by mean photon 
numbers (a^a) — \a\ 2 and (tfb) = \(3\ 2 , occupy the two 
spatial input modes of the BS. In the Fock-basis the co- 
herent state can be represented as 

|a) = £e-^" |n) = £ e -^(atr|0>, (3) 

n— n— 

and similarly for \/3). 

To account for the cases of photons being distinguish- 
able and indistinguishable at the BS we must allow for 



an additional degree of freedom in each of the spatial 
modes, e.g. polarization, frequency, or time. Thus 
we write the input state at one of the BS inputs as 
\otx, ct-i) = \o.\) ® lofa), where u\ and «2 are the coher- 
ent state amplitudes in the two orthogonal modes of the 
auxiliary degree of freedom within the same spatial mode. 
We treat the coherent state at the other BS input in a 
similar way. 

For the case in which the fields at the inputs of the BS 
are distinguishable with respect to the auxiliary degree 
of freedom, the inputs to the BS are described as being 
in the state \a, 0)|0, ft) = \a, 0) ® |0, /?), whereas in the 
case of them being indistinguishable (up to a difference 
in the mean photon number) the input fields are written 
as | a, 0)|/3, 0). 

The BS is characterized by its reflection amplitude r 
and transmission amplitude t = \f\ — |r| 2 , which cause 
the input creation operators to transform as —> t& + 
ircv and w — > ire 1 + td) . With this in hand, we can 
compute the state in the BS outputs for any combination 



of Fock states at the inputs. When the two input states 
are indistinguishable, i.e. in the same auxiliary degree of 



The above calculated output modes impinge on the sin- 
gle photon detectors (SPDs). These may be characterized 
by the probability of detecting an incident single photon. 
From this single photon detection probability rj it is also 
possible to deduce the probability of detecting a pulse 
consisting of multiple photons, keeping in mind that, irre- 
spective of the number of photons, only a single detection 
event can be generated. We write p\ (n) for the probabil- 
ity for generating one detector event given n incident pho- 
tons, and it is useful to note that it relates to the prob- 
ability po(n) of detecting nothing as pi(n) = 1 — Po(n). 
The probability for not detecting n photons is, on the 
other hand, easily computed as Po(n) = (1 — rj) n . Since 



where K\\ (n, m, j, k) should be substituted with the 
factor from Eq. Q. For distinguishable inputs we 
find a similar expression for P-q (n, to) using the factor 
Kx(n,m,j,k) from Eq. |5]). It is assumed that the de- 
tector at a given spatial output mode is equally sensitive 
to photons in both auxiliary modes, i.e. it detects the 
states \k, j) and \j, k) with equal probability. 

We are now in the position to formulate an expression 



freedom, we get [35] 



(4) 



(5) 



I 

the two detectors at the BS outputs are independent, the 
probability pn(n,m) of generating a coincidence event, 
i.e. having simultaneous detection events in each of the 
detectors, given n and m photons in one and the other 
output is simply pu{n, to) = p\ {ri)p\ (to) . Thus the prob- 
ability for a coincidence detection becomes 

Pn (n,m) = [i-(i- m) n ] [i-{i-m) m ] , (6) 

where rji and 772 are the single photon detection proba- 
bilities for detector 1 and 2, respectively. Expressing the 
coincidence detection probability in terms of Fock states 
at the BS input we have 



(7) 



I 

for the different detection probabilities given a particular 
set of coherent input fields. The probability to generate 
a detection event in both detectors, given coherent input 
fields of amplitudes a and j3, is 

oc 00 / ft om\2, 

= e e e-H 2 H^(°^prw). 

L — ' L — ' n! ml 

n—0 m=0 

(8) 



|n, 0)|m, 0) -> E E K W ( n ' m '- ? '' ^) \j + k, 0}\n + m — j — k, 0) 

j=0 fe=0 



K, ( „, m , m . e-^m-^) (») ( j Y) ( n+ ::r *)■ 

where the binomial coefficient (*) = y \^L y )\ ■ For distinguishable input fields the output state is slightly simpler 



,0)|0,m) ^EE^( n ' m '-?' 

k ) \ j, k)\n-j, m-k) 

j=0 k=Q 



K x (n,m,j,k) =J2J2 tm ~ k+j ( ir "> 

j=0 k=0 



■ ■ (AC' 1 

,-k 



Pn(n,m) = y^y^ \K\\ (n, m,j, k)\ 2 p n (j +k,n + m-j-k) 

n m 

= EEl^ii^ m ^'' fc )l 2 [i-(i-'/i) J+fe ] [i-(i-m) n+m - j - k ] 

j=0 k=0 
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(Note that to distinguish the probability in Eq. Q , which 
is applicable to Fock states, from that inEq. (pi), which 
applies to coherent state inputs, we use P to denote the 
former and V for the latter.) This allows us to derive the 
visibility of the HOM interference on the two detectors 
as 

Vn(a,/3,7?i,772,r) = — — — , (9) 

where we have spelled out the parameters that affect the 
value of the visibility. The quantity Vn is referred to as 
the HOM visibility. 



Simplified model for HOM visibility. 

To gain some intuitive understanding of the way the 
HOM visibility is affected by the experimental parame- 
ters we resort to a couple of approximations. Firstly, we 
assume equal mean photon numbers at the inputs of the 
beam-splitter, \a\ 2 — \/3\ 2 = (j,, the BS ratio to be 50:50 
(i.e. r = t = 1/V2), and the detectors to have equal sin- 
gle photon detection probability r\\ = rj 2 = tj. Secondly, 
since we normally work at very low mean photon numbers 
(j, < 1 only the first couple of terms of Eq. ^ need to be 
included. Specifically, we Taylor expand e~^/ 2 and keep 
only terms in the sum up to 2 nd order in /i. Thus, for 
the coincidence detection events we get the probabilities 



1,11 2^ 
-n_L 2 2 

which results in a HOM visibility of 
Vn = l 



(10) 
(11) 

(12) 



A key point is that the HOM visibility of 50% is inde- 
pendent of the mean photon number /i. This observa- 
tion can be explained by noting that in this low order 
treatment the coincidences in the case of indistinguish- 
able input modes stem mostly from events in which two 
photons are present at the same input, which occurs with 
probability paP2 + P2Po- For distinguishable input modes 
the coincidences stem from all events that contain two 
photons at the input, i.e. pip± +P0P2 + P2Po- Since, ac- 
cording to Eq. ([3]), for coherent input states, all of these 
probabilities scale in the same way with the mean photon 
number, their ratio, and thus the visibility of Eq. (f9l), is 
constant for all mean photon numbers. 



Compilation of experimental results for HOM 
interference at the few-photon level. 

Here we show the plots of coincidence count rates on 
which the few-photon values in Table [T] of the main text 



are based. We restate that coincidence count rates are 
proportional to coincidence probabilities by a factor that 
is given by the average experimental repetition rate. 
Moreover, when calculating the HOM visibility, only the 
relative probabilities or count rates in a measurement 
are important. In the experiments we change the mutual 
polarization, time separation, or frequency difference 
of the pulses at the BS (in the main text referred to 
as HOM-BS) input as explained in the earlier in the 
Supplementary Information. 

Deactivated memories: We present the data in order 
of the number of activated memories starting with none, 
i.e. pulses merely pass through attenuated to the BS. In 
Fig. [6^,) we show the coincidence counts as we vary the 
polarization difference of the pulses at the two inputs of 
the BS. Fitting the data to a sine function we obtain a 
visibility of V = 50.96±5.56%. In Fig.^) we display the 
coincidence counts as we step the temporal separation 
of the pulses at the two inputs of the BS. The count 
rates for these measurements are generally higher than 
all the other count rates presented. This is because this 
data was acquired by looking at coincidences between the 
transmitted part of the probe pulses in the configuration 
of two active quantum memories (shown in Fig. [8]3)). 
Hence, the balancing of the mean photon number in the 
transmitted pulses done less meticulously, which is the 
most likely reason for the observed lower visibility of V = 
42.43 ± 2.27% in this case. 

Fig. [6j;) shows the coincidence count rates as function 
of the frequency difference of the two pulses at the BS 
inputs. The horizontal line and surrounding shaded 
band shown in Fig. |6k) - as well as in Fig. R;) - give 
the coincidence counts for completely distinguishable 
input photons as obtained by making the polarizations 
orthogonal. As noted earlier in the Supplementary 
Information, it is necessary to resort to the polarization 
degree of freedom in order to make the pulses completely 
distinguishable. The visibility from the fit is noticeably 
lower than that obtained when we change the other 
degrees of freedom. There are two main reasons for 
this. The first is that, in order to generate pulses with 
different frequencies, we drive the AOM at the limits of 
its bandwidth. This, in turn, necessitates setting the 
RF drive signal amplitude high whereby the frequency 
purity of the signal is contaminated by higher-order 
harmonics. Although it is not expected to change the 
maximal interference value occurring when the pulses 
are generated with the same modulation frequency, it 
will alter the shape of the interference as a function 
of the pulse frequency difference. Hence, the fitted 
Gaussian curve, assuming a Fourier limited pulse, may 
not correctly reproduce the actual frequency dependence 
of the interference. Indeed, the minimum coincidence 
rates consistently fall below the fitted curve. A second 
factor reducing the observed visibility is related to the 
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FIG. 6: HOM interference manifested in coincidence counts between BS outputs with inactive memories, a) Changing the 
polarization angle between the pulses yields a HOM visibility of V = 50.96 ± 5.56%. b) Varying the temporal overlap of the 
pulses produces V = 42.43 ± 2.27%. c) Altering the frequency overlap of the pulse spectra results in V = 43.72 ± 1.70%. 



need to adjust the AOM drive amplitude to balance the 
bandwidth limitation. The limited accuracy with which 
we are able to estimate the appropriate RF amplitude 
results in significant scattering of the coincidence counts 
due to variations in input pulse intensities. To amend 
this we have found that it is necessary to normalize the 
points to the count rates on the individual detectors, as 
indicated on the y-axis of plot Fig. [7J;. Unfortunately, 
the manifestation of the HOM interference in the single- 
detector count rates - which will be elaborated later 
in the Supplementary Information - means that such a 
normalization procedure tends to reduce the visibility in 
the coincidence counts. 

One active memory: Next in line are the plots for the 
case in which only memory a is activated, while the 
other is left inactive. In Fig. [7] we present the coincidence 



count rates when changing the same degrees of freedom 
as in case of both memories being inactive. Addition- 
ally, in Fig. [7}i, we plot the coincidence count rates 
when changing the storage time in the quantum memory. 



Two active memories: Lastly, we present the plots for 
the case in which both memories are activated. Due to 
limitations in our current setup it is not possible to simul- 
taneously generate two quantum memories with different 
storage times, and therefore we do not acquire a storage 
time scan when both memories are active. Furthermore, 
we skip the characterization with respect to the spectral 
degree of freedom. The coincidence count data for the 
remaining two degrees of freedom are plotted in Fig. [8j 
which also includes the appropriate fits. 
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FIG. 7: HOM interference manifested in coincidence counts between BS outputs with one active memory, a) Changing the 
polarization angle between the pulses yields a HOM visibility of V = 55.51 ± 4.09%. b) Varying the temporal overlap of the 
pulses produces V = 47.57 ± 2.96%. c) Altering the frequency overlap of the pulse spectra results in V = 42.40 ± 3.51%. d) 
Varying the storage time of the quantum memory and thus the temporal overlap of the pulses yields V = 44.4 ± 6.9%. 



Manifestation of HOM interference in single 
detector counts. 



We also evaluate the effect of the two-photon interfer- 
ence on the counts registered by a single detector. This 
is easily done by amending the detection probability to 
the case of one detection event in one detector and any 
number of events x (i.e. x — 0, 1) in the other detector. 
We arrive at 



p lx (n,m) = 1 - (1 - r)xY 



(13) 



This expression can be inserted into Eq. ^ to cal- 
culate Pj'j^^n, m), which, through Eq. (JsJ) , gives us 
i a i P)' ana - horn which the single- detector visibil- 
ity Vix is defined analogous to Eq. Efy. 



We can formulate a simplified expression by using the 
same approximations as in the case of coincidence detec- 
tions: 



(14) 
(15) 

(16) 



from which we get the single-detector visibility 

V = W 
lx 4 + 2(4-?7) M ' 

In the limit of low detector efficiency, Vi x ~ 0, since, in 
that case, the probability of detecting two photons im- 
pinging on the detector is simply twice that of detecting 
one. This nulls the limitation that only a single detec- 
tion event can be generated per pulse. Furthermore, the 
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FIG. 8: HOM interference manifested in coincidence counts between BS outputs with two active memories, a) Changing the 
polarization angle between the pulses yields a HOM visibility of V = (53.1 ± 5.3)%. b) Varying the temporal overlap of the 
pulses produces V = (46.1 ± 3.2)%. 



single-detector visibility also goes to zero for very low 
mean photon numbers. In this case it is very unlikely 
to have two photons either at the same or at different 
input ports of the BS, hence most of the single detector 
counts stem from single photons from either one or the 
other input of the BS. It is interesting to note that if rj 
is known for a detector, then, from observing the single- 



Experimental results on HOM interference 
manifested in single-detector counts 



detector visibility (see Eq. ( 16 1 ), it is in principle possible 



to estimate the mean photon number per pulse //. 
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FIG. 9: Plots of single-detector visibility as a function of the 
mean photon number for detectors with a range of single pho- 
ton detection probabilities n. The n — 70% trace, highlighted 
with a dashed line, corresponds approximately to our detec- 
tors, which have 65% < n < 75%. 

Another important consequence of the manifestation 
of two-photon interference in the single-detector counts is 
that the single-detector counts cannot generally be used 
to normalize the coincidence counts w.r.t. fluctuations in 
the input pulse intensities. Only for detectors with low 
detection efficiency or very low mean photon numbers, in 
which case Vi x ~ 0, is this normalization possible. 



First, in Figure 10 



we present the single-detector 
counts corresponding to the coincidence counts depicted 
in Figure [7^i,b. In the case where we vary the polariza- 
tion and time separation we see a clear change in the 
single-detector counts, which, moreover, is evidently cor- 
related with the change in coincidence counts. The count 
variation due to the two-photon interference is somewhat 
masked by the single-detector count scatter, which is due 
to intensity fluctuations mainly in the light going through 
the 10 km delay line. We fit the data in Figures [TOfc, and 
b with a sine and Gaussian function, respectively. For 
the former we find a mean photon number of /1 = 0.52 
while from the latter we estimate \x = 0.54. From the 
number of single-detector counts there is some evidence 
to conclude that the light intensity is about 15% higher. 
To this should be added about 25% uncertainty for the 
intensity at the BS w.r.t. the intensity at the detector 
due to variation in the loss in the fibre mating sleeves. 
Finally, the scatter of the counts makes the fits them- 
selves rather uncertain. Nevertheless, the mere fact that 
the two-photon interference is manifested in the single- 
detector counts validates the order of magnitude of the 
mean photon number, as depicted in Fig. [9j 

Figure 11 depicts the single-detector counts corre- 
sponding to the coincidence counts depicted in Fig- 
ure[8^i,b. Again, from fitting the appropriate functions to 
the polarization and time data yields visibilities around 
7%, corresponding to mean photon numbers of around 
fj, = 0.5. 



13 




40 80 120 160 200 240 -30 -20 -10 10 20 30 

polarization angle [deg] pulse time separation [ns] 



FIG. 10: HOM interference manifested in single-detector counts in the case of one active quantum memory when changing a) 
polarization and b) time difference between pulses at BS input. For the polarization scan in a) we find Vi x = (7.51 ± 3.80)% 
and for the time scan in b) we get Vi x = (7.75 ± 3.25)%. For this measurement we only recorded the single-detector counts 
from Si-APD 1. 
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FIG. 11: HOM interference manifested in single-detector counts in the case of two active quantum memories when changing a) 
polarization and b) time difference between pulses at BS input. For the polarization scan in a) we find Vi x = (8.64 ± 2.50)% 



and Vi x = (7.60 ± 2.36)% for Si-APD 1 and 2, respectively. 
Vix = (6.23 ± 1.61)% for Si-APD 1 and 2, respectively. 



For the time scan in b) we measure Vi x = (6.38 ± 2.01)% and 



Bell-state measurement. 

In this section we derive an analytical expression for 
the coincidence count rates corresponding to projections 
onto the \ip~) Bell-state for time-bin qubits detected by 
the two detectors at the output of the HOM-BS. To that 



end, we will introduce a number of approximations as we 
did previously in order to calculate the HOM interfer- 
ence in the coincidence counts. In the limit of low mean 
photon numbers, two coherent states impinging onto the 
two inputs of a 50:50 BS can be represented in terms of 
Fock states as 
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|V>>«6 = VpO-, l)|H)a,6 + v / P(2,0)|20) a , 6 + v'p(0,2)|02) o , 6 

1 



(17) 



V2! 



V^O)((at) 2 ® /) + Vp(P^)(I ® (S f ) 5 
I 



|00) O)i 



where the subscripts on the state vector refer to the or- 
der of listing the input modes, i.e. |00) n ,b = |0}„ <g> \0)b- 
The factors written as p(n, m) denote the probability 
of having n and m photons in mode a and b, and 
are given by p(n,m) = \( a (n\ ® b (m|)(|a) a ® \P)b)\ 2 = 

Stemming from the low mean 



(M 2 ) 



2\nl 



\2\rn 



photon number assumption, we do not include terms with 
more than two photons. Assuming that our detectors are 



noiseless, terms with a total of one or no photons are also 
left out as they cannot generate any coincidence counts. 

For a time-bin qubit, the Fock state is created in a 
superposition of two temporal modes, i.e., an early (e) 
and a late (I) mode, by the creation operators for the 
spatial input mode (x^ — a^.tf) of the beam-splitter, 

as 



|00) a 



(18) 



where cos(^) and sin(^) are the amplitudes of, and </> x 
is the relative phase between, the two temporal modes 
composing the time-bin qubit. The subsript xe refers to 
the early time-bin of the spatial mode x and similarly 
for xl. Note, that we sometimes simplify the notation 



for the time-bin qubit states as \e) x = 1 10) 



xe,xl 



= r*t 



(x\® 



If we insert the expression in Eq. ( 18 ) in 



I)\00) xe 

place of the a and b operators in Eq. (1171) we get the 



expression for the wavefunction \ip{0 a , <f> a , 9 b , 4>b))ab f° r 
time-bin qubits at the HOM-BS inputs. We split this 



expression into the various contributions given in Eq. ( 17 ) 



((tf) 2 ®I)\00) ab ^ 



* cos(|) sin(|) + ie«- sin(|) cos(f ) J + did] 

-cos(|)sin(|)-e-sin(|)cos(|))(^_ a t d]) 

+ ie ^+« sin (|) sin (|) ((4f + (di) 

+ *cos(|)cos(|) ((ct) 2 + (4) ; 



|0000) 



ce,cl,de,dl 



2e*-co B (£)sin(£) (ctct-dldj 
^ 2e ^ cos (|) sin (|) (cid}+c}dl) 
(|) ((ct) 2 + i244-(4) 



cos 



{c\f + i2c}d]-{d}) 



|0000) 



ce.cl.de.dl 



(19a) 



(19b) 



and similarly for (I <£> (6^) 2 )|00) a f,. Again, the subscripts poral and spatial modes, e.g. ce labels the early bin of 
on the state vector refer to the order of listing the tern- 
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the spatial output mode c. 

We will look for coincidence detection events that 
correspond to projections onto the Bell-state |'0-)cd = 
- c\dl)\WQQ) cetC i.de,di- Such projections corre- 
spond to a detection event in the early time-bin in one 
detector followed by a detection event in the late time-bin 
in the other detector. This projection occurs with a prob- 



ability V-(0 a ,(/)a,0b,(l>b) = | cd{lp-\^{O a ,(j) a ,db,(j>b))cd\ 



which can be computed by combining Eq. ( 19 ) with 
Eq. (17). Assuming equal mean photon numbers at the 
2 = |/3 1 2 = n and averaging over the coher- 
i.e. the complex angle between a and 



two inputs | a 
ent state phases 
/3, we get the expression 



V-{9 a ,4> ai 6 b ,4> b ) oc 



^ 2 e~ 2 ^ 



4 sin 



2 / "« 



+ sin 2 (0 a ) +sin 2 (0 b ) 



- 2sin(6» a ) sin(0 6 ) 1 + cos(^ a - 



(20) 



With this we are able to calculate the probabilities of 
projection onto for different combinations of qubits 
at the two BS inputs, i.e. for different choices of the 
angles 9 X and <p x . In turn, this allows us to calculate the 
Bell-state measurement error rate as 



(21) 



where V_ is the projection probability when the two 
input qubit states are identical, i.e. <j) a — (f>b and 9 a = 9b, 



while V^z is the projection probability for two orthogonal 
input qubit states. This is also defined in terms of count 
rates in Eq. ^ in the main text. We will now treat a 
number of relevant cases. 



Expected and observed error rates when 



= 0. 



Using the simplified notation this corresponds to the 
case were the input qubit states are of the form = 
cos(^)|e) +sin(%-)|/). When depicted on the Bloch 
sphere these qubits span the xz-plane. Using Eq. (20) 
we compute the projection probability as 



7>_(0 Q ,OA,O)a 



M 2 e- 2 ^ 



4 sin 



?a + 0b 



sm 2 (6 a ) + sin 2 (0 fc ) - 4sin(0 a ) sm(9 b ) 



(22) 



We are interested in the probability V_ for the case in 
which the input qubits are parallel (9 a — 9b) and Vzz for 
the case in which the input qubit states are orthogonal 
(6 a = 9j, — 7r). Specifically, when we prepare two qubits 
(one at each input of the BS) in state |e), or two qubits in 
state |^), we expect V_ = 0. The probability for observ- 
ing a projection onto increases as we change 9 a (or 
9},), and reaches a maximum T>z\ if one qubit is in state 
|e) and the other one in \l). Hence, using the expression 



for the error rate above 



(Eq. (2lJ), 



we find e 



(att) _ 
e/l 



0. 



We now turn to measuring the coincidence rates for 
all combinations of |e) and \l) input states, and thus 
extracting v\_ and V^, using 0.6 photons per qubit 
at the memory input. More precisely, we prepare the 



input qubit state 



|e)b to measure V_ 



IKD 



and 



then \e) a <£) \l)b to measure V. 



-L(l) 



Subsequently, we 



prepare the input qubit state |^) a ® \l)b to measure 



and then |^) a (g) \e)b to measure V±^> . These 



,-L(2) 
l(2h 



yield the average values Vl = (V Ul) +Vl (zl )/2 and 
= [V^ 1 ^ + V^ 2 ^) /2, from which we compute the 
experimental error rate e^ p ^ = 0.039 ± 0.037, which is 



"e/l 

near the theoretical lowest value of e 



(att) 



0. 



Expected and observed error rates when 9 a = 9b = n/2. 
In this case the two input qubits are in equal super- 
positions of early and late bins, that is of the form 
|^) = ^(|e) + e^\l)). On the Bloch sphere these are 
qubits that lie in the xy-plane. In this case we compute 

P_(7r/2,^,,7r/2,&) cx (2 - cos(0 a - fa)) , 

(23) 

Thus the \ip~) Bell-state projection probability is small- 
est - but nonzero - when <j) a — 4>b — 0, i.e. the qubit 
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states are parallel, and largest when the phases differ by 
7T, i.e. the qubit states are orthogonal. Inserting these 
values for v\_ and Vz: into Eq. ( 21 ) results in an expected 
error rate of e^ 1 '^ = 0.25. 

Using again 0.6 photons per qubit, we measure 
the coincidence counts for <j> a — fa = and 7r giving 
us p! and P^z , respectively. From these we get an 
error rate of = 0.287 ± 0.020, which is slightly 

above the theoretical bound. This indicates that either 
the measurement suffers from imperfections such as 
detector noise or the modes at the BS are not com- 
pletely indistinguishable, which in turn could be due 
imperfectly generated qubit states or imperfect storage 
of the qubit in the quantum memory. To be conser- 
vative in our assessment of our quantum memory we 
assume that the entire increase of the measured values 
of e( oxp ) is due to the memory fidelity being less than one. 

Bounds for attenuated laser pulses stored in quantum 
and classical memories: We now compare the perfor- 
mance of our Bell-state measurement to a number of 
relevant bounds assuming always that any imperfections 
arise from the imperfect storage of the photon in the 
memory. We will derive bounds to the error rate in the 
case of one qubit being stored in either a classical mem- 
ory (CM) or quantum memory (QM). To accommodate 
this scenario we assume that the memory performs the 
following operation |^)(^| -> F\^) (ip\ + (1- F)^) (i/j^l, 
where F denotes the fidelity of the stored state and \ip' L ) 
is the state orthogonal to \if>). For a classical memory 
F CM = 0.667 [47] whereas for a quantum memory 

F QM = 1 

Doing the replacement v\_ Fvl + (1 - F)Vz: and 
likewise for Vzz we can express the error rate expected 
after imperfect storage of one of the pulses partaking in 
the Bell-state measurement: 



bounds e 



(att.QM) 



— and e ^* t ' CM ) = 0.333. This clearly 
shows that a classical memory would cause a larger er- 



FVl + (1 - F)V± 



(24) 



where in this case the probabilities V_ and Vzz refer to 
those expected without the memory. Since the expected 
values for "Pi and Vzz differ between the e/l and +/— 
bases we treat them separately. 



Beginning with the e/l basis we use Eq. (24) with the 
values from Eq. ( 22 ) to derive a bound for the error rate 



of the Bell-state measurement for one of the two qubits 
being recalled from a quantum or a classical memory. We 



find that e 



(att) 

■■II 



ror rate than the e 



(exp) 
:/l 



0.039 ± 0.037 measured after 



storage in our memory. We can also reverse the equa- 
tions and estimate our memory's fidelity based on the 



(exp) 
e/l 

Eq. p4|, we deduce the value F^jf = 0.961 ± 0.037 



measured error rate. In this case, inserting e 



into 



We now turn to the +/— basis. For attenuated laser 
pulses we insert into Eq. (24) the values V_ = 1/4 
and Vz: = 3/4 computed from Eq. (23), which enables 



us to relate the error rate to the memory fidelity as 
_ = (3 — 2F)/4. Thus, one obtains the theoretical 



lower bound on the error rate e 



(att.QM) 
+/" 



0.250 for an 



0.417 



ideal quantum memory (_pQ M = 1) and g^**^ 1 ^) 
with an optimal classical storage device (P CM = 2/3). 
We make the observation that our experimental error 
rate = 0.287 ±0.020 is much below the bound for a 

classical memory. Based on the experimental error rate 
= 0.287 ± 0.020 we derive an experimental value 
for the memory fidelity of F™P = 0.926 ± 0.041. The 
estimates of the memory fidelity F°J^ and F^J^ derived 
from our measurements in two bases are equal to within 
the experimental error. This together with the fact that 
their values are well above 0.667 reaffirms our claim that 
our storage device outperforms a classical memory. 

We emphasize once more that we have assumed that 
the reduction in error rates is due solely to the memory 
and thus indicates the fidelity of the memory. However, 
this is likely not the case as imperfections in the state 
preparation and detector noise also contribute to the 
reduction in error rate. 

Bounds for single photons stored in quantum and clas- 
sical memories: Although we do not use single photon 
sources for the experiments reported here, it is inter- 
esting to determine how well our results measure up to 
those that could have been obtained if single photon 
sources had been employed. In the following we will 
derive the error rate for the Bell-state measurement 
using qubits encoded into single photons. To this end we 



step back to Eq. (17), and note that for single photon 
sources all probabilities are except for p(l, 1), which 
describes the probability of having a single photon at 
each BS input. Thus, in the output state we only need 
to keep the terms from Eq. (19a), which in turn means 



that the Bell-state projection probability can be written 



1 — F, and hence we establish the two as 



P-(0a,4>a,0b,4>b) oc 



1 



sm 



h + o b 



sm 



- sin(0 a ) sin(6> fc ) cos(fa- fa) 



(25) 
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It is easily seen that for any two parallel input qubit 
states {9 a = &b and 4> a — 4>b) we get pI = 0. There- 
fore, irrespective of the projection probability for orthog- 
onal input qubit states the expected error rate is always 
e (smg) _ w } icrc sing identifies this value as belonging 
to the single photon case. 

Gauging the effect of storing one of the single pho- 
tons partaking in the Bell-state measurement in a 
mem ory is thus independent of the basis and using 



Eq. [ 

P (sing,i 



we derive e ( sin s<Q M ) 



1 



F QM 



and 



= 0.333. Contrasting the error rate expected 
for a photon stored in a classical memory with the two 



values 



,( cx p) 

'e/l 



0.039 ± 0.037 and e 



(exp) _ 
+ /- 



0.287 ±0.020 



obtained experimentally, we recognize that both are well 
below e^ slng,CM ' > . This means that even with a single 
photon source at ones disposal the error rates that we 
measured could not have been attained with a classical 
memory. 



visibility as 



V 



px _p\\ 



(26) 



analogous to Eq. (JlJ in the main text. Using val- 
ues obtained from a cosine fit to the data in Fig. [T2| 
yeilds V+T = (62.9 ± 5.2)%. Comparing Eq. (l26|) with 
Eq. (21| it is easily seen that V and e are related as 



e = (1— V_|_/_)/(2— V+/_). We can then use the expected 
error rates to find the corresponding Bell-state measure- 
ment visibilities. Using e°**_ = 0.25 we get a theoretical 
value V^y_ = 66.7%. In conclusion, our experimental 
Bell-state measurement visibility is only slightly below 
and within the experimental error actually equal to the 
expected value. 
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FIG. 12: Rate of projection of pairs of time-bin qubits with 
relative phase <j> a — (j)^ onto \ip~)- Each data point was ac- 
quired over 60 s 



incidence detection probability reaches its maximum V_ 
when two input qubits are orthogonal (<p a — §b = tt) and 
when they are identical (0 a — = 0) it reaches a mini- 
mum V_ ■ It is natural to define a Bell-state measurement 
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